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Abstract. - In order to define the structure of trophic network of Ichkeul Lake in Tunisia, diet of main teleost 
fish species was determined during two seasons, wet and dry. A total of 491 stomachs from 16 teleost species was 
analyzed. The analysis revealed spatial variability (East under marine influence and West with more continental 
effect) and temporal one (wet and dry seasons). Based on stomach contents, nine trophic groups were revealed, 
including eight monospecific groups. The main prey are mudflat snails ( Hydrobia ventrosa, H. acuta), bivalves 
(Abra alba, Cerastoderma glaucum), amphipods and isopods, as well as seagrass ( Potamogetonpectinatus, Rttp- 
pia cirrhosa ) and water brackish algae (Ulva spp. and Chaetomorpha spp., Cladophora sp., Polysiphonia sp.). 


Resume. - Etude du regime alimentaire de l’ichtyofaune du lac Ichkeul (Tunisie). 
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Afin d’etablir la structure du reseau trophique du lac Ichkeul en Tunisie, le regime alimentaire des principales 
especes de poissons teleosteens pendant deux saisons, humide et seche, a ete determine. 491 individus apparte- 
nant a 16 especes de poissons ont ete analyses. L’analyse des contenus stomacaux arevele une variabilite spatiale 
(Est sous influence marine et Ouest plus continentale) et temporelle (saison humide et saison seche). Sur la base 
des contenus stomacaux, neuf groupes trophiques ont ete crees, dont huit sont mono-specifiques. Les principales 
proies sont les gasteropodes des vasieres ( Hydrobia ventrosa, H. acuta), les bivalves ( Abra alba, Cerastoderma 
glaucum), les amphipodes et isopodes, ainsi que les phanerogames ( Potamogeton pectinatus, Ruppia cirrhosa) 
et les algues des eaux saumatres ( Ulva spp. et Chaetomorpha spp., Cladophora sp., Polysiphonia sp.) 


The trophic structure of ichthyofauna communities is dif¬ 
ficult to characterize and to predict because of the diet plas¬ 
ticity according to the availability of food resources (Piet, 

1998) . The partition of trophic niche was recognized as head 
main factor in the communities structuring (Schoener, 1974; 
Ross, 1986; Sibbing et al., 1998). Fish diet informs on the 
presence, the abundance and the availability of prey, and 
also on trophic potential of the ecosystem. The strategies 
of searching food and variations of the species diet from an 
environment to another and/or from a season to another one, 
inform about the fish adaptation capacities in front of the 
environment constraints (Paugy and Leveque, 1999). The 
analysis of stomach contents allows taxonomic identifica¬ 
tion of the various groups and species belonging to the struc¬ 
ture of trophic network and the trophic interactions within 
the complex ecosystems (Piet, 1998; Paugy and Leveque, 

1999) . 

The importance of the study site lies in its specificity; the 
Ichkeul Lake represents one of the three main components 
of the Ichkeul Park. It is a protected ecosystem recognized 
for the value of its biological heritage, in particular birds and 
fish species (BCEOM, 1995). 


Several studies on the lacutrine Ichkeul ecosystem have 
focused on specific richness and bioecology of some fauna 
or floral species (BCEOM, 1995; Chaouachi and Ben Has- 
sine, 1998; Kraiem et al., 2003; Casagranda et al., 2005b; 
Sellami et al., 2010) but none was interested in the structure 
and trophic functioning of this ecosystem (except for two 
works by Casagranda et al. (2005b) and Casagranda and 
Boudouresque (2010) that did present a simplified concep¬ 
tion of the general functional model of Ichkeul ecosystem). 
The knowledge of ingested preys and feeding habits of fishes 
in any environment is essential to understand their adapta¬ 
tion to the environmental drivers. This work aims to study 
spatiotemporal variability of ichthyic community diet of the 
Ichkeul Lake. 


MATERIALS AND METHODS 

Study site 

The Ichkeul Lake (80 km 1 2 ) is situated in the north of 
Tunisia (37°8r00”N; 9°83’30”E). Its area varies from a 
season to another, from 12,000 ha in winter to 8,000 ha in 
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summer (ANPE, 2007). It is supplied by fresh water from 
a surrounding watershed of 8,000 km 2 through seven oueds 
(rivers), and it is connected by Tinja canal along 5 km to 
the Bizerta lagoon, which opens on the Mediterranean Sea 
(Fig. 1). The average depth is 1 m with a maximum of 2 m 
registered in winter period. 

In summer, the high evaporation lowers the water level, 
which allows water of Bizerta lagoon dumping in the lake. 
The salinity shows considerable seasonal variability from 3 
near the oued outfalls during springtime, to more than 45 on 
the Tinja channel (canal) during autumn season. 

Salinity and temperature were measured on the site dur¬ 
ing sampling trips using portable multimeter (Thermo Orion 
meter, model 115A+). 

The sampling was designed as to cover whole lake sec¬ 
tors, by considering the oued embouchures, the variation of 
water level, the influence of fresh and/or marine water and 
finally the variability of macrophytes cover. Samples were 
performed in two seasons from March 2011 to April 2012; 
dry season from June 2011 to November 2011, wet season 
from March 2011 to May 2011 and from December 2011 
to April 2012. For stomach content analysis, the fish speci¬ 
mens were collected by means of benthic dragnet (meshes of 
10 mm) and trammel nets (meshes of 26/28 mm and 30/40 
mm). These nets are deployed in the morning and recovered 
at the end of the day. 

Additional commercial fish samples were collected 
from catches of “Tunisie lagunes company”, the only owner 
exploiting fishes of the lake. The fishing means vary between 
trammel nets (mesh of 30/40 mm), traps and fyke nets 
(capechades) (mesh of 6 mm), fixed gillnet (eels dams) and 
fixed fisheries (weir deployed nearing Tinja canal sluice). 


Nets are deployed in the afternoon before sunset and 
are recovered the following day in the morning. Nets are 
deployed monthly to have seasonal samples covering all 
over the lake surface, as much as possible. The wet season 
is considered from December to May, while the dry season is 
defined from June to November. 

The small sizes species are sampled by means of dragnet, 
seine net and larvae net. As possible, from 15 to 20 individu¬ 
als per species are collected to obtain 10 full stomachs. 

Sample processing 

The sampled individuals were beforehand weighed and 
their standard length (FS) was measured to the nearest mil¬ 
limetre. In the laboratory, each stomach was opened into a 
Petri dish before being kept preserved in ethanol solution. 
The observation of the stomach contents was realized under 
a binocular microscope. Preys were sorted and identified to 
the species level when possible. 

Stomach content analysis 

Several methods can be used for fish stomach content 
analysis. The choice of these methods must be in agreement 
with the studied objective and the type of required data (Hys- 
lop, 1980). 

For our study, the volumetric method of estimation was 
used (Hynes, 1950; Rounsefell and Everhart, 1953; FAO, 
1974). Within the study framework, we opted for a volumet¬ 
ric analysis method, the method of points (Rounsefell and 
Everhart, 1953; FAO, 1974). It is a variant of the visual esti¬ 
mation method. Instead of directly evaluating volume using 
eyesight, a number of points were given for each item (prey) 
found in the stomach according to its volume. The category 
having the highest volume in the stomach content will have a 


Figure 1. - Geographic location of 
Ichkeul Lake: its water network and its 
communication with Bizerta lagoon. 
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maximum of 16 points. Other items will have points varying 
between 8,4,2,1 and 0 according to their volume compared 
to the major prey category. 

The volume percentages of each prey were then calcu¬ 
lated (Hynes, 1950; Rounsefell and Everhart, 1953; FAO, 
1974). For each prey item, its volume percentage (%V) was 
calculated using the following formula: % V prey volume 
(i) = [number of points for prey (i)/(number of points for all 
preys)] xlOO. 

The results of the stomach content analysis were quan¬ 
tified using numerical indexes (Hynes, 1950; FAO, 1974; 
Hyslop, 1980). The most known index is the Frequency 
Index (IF): 

IF = Ni/N x 100, where Ni is the number of stomach contents 
for the species i, and N is the number of non-empty analyzed 
stomachs. 

For each fish species, we expressed the volume percent¬ 
ages (%V) of preys, which constitute its diet. 

The vacuity coefficient (Cv) is the ratio of the number of 
empty analyzed stomachs N e relative to the total number of 
analyzed stomachs by each species N t . Cv = N e /N t , where N t 
is the total number of analyzed stomachs per species and N e 
the number of empty analyzed stomachs. 


Statistical analysis 

Diet variations were studied according to the sector and 
season sampled for each species. 

Diet similarities between the different fish species were 
evaluated using the simplified Morisita index, proposed 
by Horn (1966), which is also known as the Horn-Morisita 
index (Krebs, 1999). This similarity index was calculated 
using the reconstruction of the weight data using the follow¬ 
ing formula: 


Cmh 


2 £f=i Pa,i x Pb,i 

Ef=i Pa/ + £f=i Pb,i‘ 


C,„h is the similarity index of Horn-Morisita between the 
predator species A and B; S is the total number of prey spe¬ 
cies identified in the two predators diets, P A l is the reconsti¬ 
tuted average weight proportion of prey i consumed by spe¬ 
cies A and Pbj is the reconstituted average weight proportion 
of prey i consumed by species B (Potier et al., 2007). The 
similarity index varies between 0 (no preys in common) and 
1 (complete similarity) for a similarity threshold generally 
set at 0.6 for a “significant” similarity between two species 
(Zaret and Rand, 1971). C m /, was calculated within a 95% as 
confidence interval. 

The diet similarity matrix of Cindex was completed 
by a Hierarchical Ascendant Classification (HAC) using 
1 - C m h (Euclidian distance, single link) in order to identify 
groups having similar diets (Zaret and Rand, 1971). 


A Correspondence Analysis (CA) (Benzecri, 1973; 
Benzecri and Benzecri, 1984) based on the abundance and 
occurrence of prey species was used to examine the trophic 
tendencies of the sampled species. The aim of this analysis 
was to find groups with similar diet affinities. The CA, the 
similarity indices and the HAC were performed using Past 
software version 3.0.1 (Hammer et al., 2001). 


RESULTS 

Environmental characteristics 

The annual temperature of Ichkeul Lake varies between 
33.8°C in summer and a minimum of 10.1°C in winter. 
There is some thermal variation between eastern and western 
sectors of the lake. The temperature is slightly hotter in the 
eastern sector during the months of July, August and Sep¬ 
tember, as well as months of February and April (Fig. 2A). 
The annual salinity oscillates between salty marine waters 
during hot season and freshwater during cold season, with 
autumn maximum at 37.5 and winter minimum at 0.2. The 
eastern and western parts are characterized by highly differ¬ 
ent salinity from May to September, with eastern zone hav¬ 
ing the highest salty rates (Fig. 2B). 




Figure 2. - Evolution of the average temperature (°C) (A) and salin¬ 
ity (B) in both East and West sectors of Ichkeul Lake during the 
sampling period. 
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Stomach content analysis 

A total of 491 individuals belonging to 16 Actinopterygii 
fish species were analyzed during this study. Identification 
has allowed us to recognize nine marine species (i.e. having 
a part of life cycle in the sea; Anguilla anguilla (Linnaeus, 
1758) (Ana), Dicentrarchus labrax (Linnaeus, 1758) (Dil), 
Belone belone (Linnaeus, 1761) (Beb), Solea senegalensis 
(Kaup, 1858) (Sos), Engraulis encrasicolus (Linnaeus, 1758) 
(Ene), Mugil cephalus (Linnaeus, 1758) (Muc), Liza rama- 
da (Risso, 1827) (Lir), Liza aurata (Risso, 1810) (Lia) and 
Liza saliens (Risso, 1810) (Lis), six sedentary species (i.e. 
living their entire life cycle within the lake; Aphanius fas- 
ciatus (Valenciennes, 1821) (Apf), Atherina boyeri (Risso, 
1810) (Atb), Pomatochistus microps (Krpyer, 1838) (Pom), 
Syngnathus abaster (Risso, 1827) (Syab), Syngnathus acus 
(Linnaeus, 1758) (Syac) and Syngnathus typhle (Linnaeus, 
1758) (Syty) and one freshwater species (i.e. living a part of 
their life cycle in the rivers; Barbus callensis (Valenciennes, 
1842) (Bac)) were represented in our sampling. 


The 16 species were collected during the two seasons, but 
in different quantities. 256 and 235 stomachs were analyzed, 
during the dry season since June 2011 to November 2011, 
and during the wet season since December 2011 to May 
2012, respectively. The annual vacuity rate Cv for all species 
varied between 0% (L. ramada, E. encrasicolus, P. microps, 
A.fasciatus) to 33% (A. anguilla). This index showed a high 
variability between species, sectors and seasons (Appendix 

I). 

Fish diets 

The Garfish needle ( B. belone ) was characterized by 
diet based on insects (%V > 30%) and peracarid crustaceans 
(%V = 40%). The diet of D. labrax was more diversified 
with the dominance of animal fraction essentially consti¬ 
tuted by fish (%V = 30%) and animal debris (%V > 25%). 
The anchovies were characterized by a zooplankton domi¬ 
nated diet (%V > 70%). P. microps presented a more diverse 
diet composed basically by zooplankton (%V = 45%), zoo¬ 
benthos (%V = 23%) and amphipod/isopod crustaceans 
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Figure 3. - General annual evolution in the volumetric percentages (%V) of the prey categories found in the stomach contents of Ichkeul 
fish species (the number of non-empty stomachs analyzed for each species is indicated above each species. Ana: Anguilla anguilla, Apf: 
Aphanius fasciatus, Atb: Atherina boyeri, Bac: Barbus callensis, Beb: Belone belone, Dil: Dicentrarchus labrax, Ene: Engraulis encrasi¬ 
colus, Lia: Liza aurata, Lir: Liza ramada, Lis: Liza saliens, Muc: Mugil cephalus, Pom: Pomatochistus microps, Syab: Syngnathus abas¬ 
ter, Syac: Syngnathus acus, Syty: Syngnathus typhle and Sos: Solea senegalensis. 
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Table I. - Horn-Morisita index of diets similarity between fish species sampled in Ichkeul Lake. Signifi¬ 
cance reported to C,„i, >0.6. 



Ene 

Pom 

Syt 

Beb 

Dil 

Sos 

Bac 

Apf 

Atb 

Muc 

Lir 

Lis 

Lia 

Ana 

Syab 

Syac 

Ene 

1 
















Pom 

0.84 

1 















Syt 

0.69 

0.84 

1 














Beb 

0.02 

0.34 

0.55 

1 













Dil 

0.14 

0.17 

0.22 

0.45 

1 












Sos 

0.17 

0.39 

0.63 

0.63 

0.47 

1 











Bac 

0.16 

0.32 

0.49 

0.55 

0.51 

0.91 

1 










Apf 

0.29 

0.53 

0.67 

0.53 

0.52 

0.85 

0.71 

1 









Atb 

0.11 

0.35 

0.67 

0.62 

0.41 

0.86 

0.72 

0.84 

1 








Muc 

0.29 

0.37 

0.33 

0.15 

0.40 

0.61 

0.75 

0.48 

0.44 

1 







Lir 

0.32 

0.35 

0.41 

0.19 

0.47 

0.72 

0.85 

0.58 

0.60 

0.89 

1 






Lis 

0.26 

0.31 

0.38 

0.20 

0.40 

0.70 

0.86 

0.53 

0.59 

0.92 

0.98 

1 





Lia 

0.23 

0.22 

0.36 

0.21 

0.41 

0.74 

0.87 

0.52 

0.58 

0.80 

0.96 

0.95 

1 




Ana 

0.12 

0.09 

0.29 

0.27 

0.40 

0.61 

0.78 

0.40 

0.47 

0.52 

0.72 

0.71 

0.80 

1 



Syab 

0.45 

0.46 

0.53 

0.21 

0.62 

0.62 

0.62 

0.77 

0.67 

0.58 

0.77 

0.68 

0.68 

0.51 

1 


Syac 

0.75 

0.62 

0.69 

0.14 

0.31 

0.50 

0.56 

0.47 

0.43 

0.54 

0.74 

0.66 

0.73 

0.58 

0.73 

1 


(%V > 15%). S. typhle showed a diet essentially composed 
by zooplankton (%V = 38%) and amphipod/isopod crusta¬ 
ceans (%V > 35%) (Fig.3). S. acus fed on almost equal quan¬ 
tity on zooplankton (%V = 40%) and vegetals (macrophytes/ 
algae; %V = 42%), whereas S. ahaster ingested a more 
important animal fraction (%V > 70%) mainly composed by 
animal debris (%V > 30%), zooplankton (%V = 20%) and 
micro-gastropods (%V = 10%) (Fig. 3). 

The food diet of A. anguilla was mainly composed by 
bivalves (%V = 30%), macrophytes/algae (%V > 35%) 
and to a less extent by crustaceans (%V = 10%) and fish 
(%V = 10%). Mullets were characterized by a common 
trophic spectrum for the four Mugilidae species (L. aurata, 
L. ramada, L. saliens and M. cephalus). This food spec¬ 
trum was essentially composed by macrophytes/algae 
(%V = 20-40%), zoobenthos (%V = 10-15%) and gastropods 
(%V > 10%). The ingested sediment fraction showed differ¬ 
ent contribution in each species and varied between 10% ( L. 
aurata ) to 25% (M. cephalus). 

The Barb B. callensis showed a diet which mainly includ¬ 
ed peracarid crustaceans (%V = 15-35%), various animal 
debris (%V = 6-20%) and a vegetal fraction (%V = 10-25%). 
The Big-scale sand smelt A. boyeri exhibited a prey gastro¬ 
pods contribution ( Hydrobia spp.) of 20%, while the South 
European Toothcarp A. fasciatus was characterized by a 
diet with less gastropods (1%) and amphipod/isopod crus¬ 
taceans (%V = 25%) contributions (Fig. 3). The S. senega- 
lensis diet was the more diverse, with an important predation 
on benthic crustaceans (%V = 15-40%), Nereis diversicolor 
(%V = 7%), insects, hydrobies and meiofauna. 

The different sampled and identified preys in the fish 
fauna stomach contents are detailed in table II. Their volu¬ 


metric prey contributions (%V) and occurrence coefficients 
(IF) are summarized in table III. 

Trophic groups 
Inter-species diet similarities 

At the threshold of 0.60 for Horn-Morisita index, the 
HAC revealed nine trophic groups among them, eight 
were monospecific: B. belone, E. encrasicolus, P. microps, 
S. typhle, S. acus, S. abaster, A. anguilla and D. labrax 
(Fig. 4). B. belone had a carnivorous diet with an affinity for 
insects, while E. encrasicolus was zooplanktonophagous and 
P. microps was carnivorous with an affinity for zooplankton. 

The trophic categories resulting from this stomach con¬ 
tent study (Figs 5, 6) were different from those mentioned 
in the other studies (Appendix II). The environment, trophic 
potentialities, and availability of preys influence the trophic 
behaviour of fish species in the living space. B. belone had a 
carnivorous diet with an affinity for insects, while E. encra¬ 
sicolus was zooplanktonophagous and P. microps was car¬ 
nivorous with an affinity for zooplankton. 

Among Syngnathidae, the diet of S. typhle was very dif¬ 
ferent from the other representative of its family, with an 
affinity for amphipods/isopods/zooplankton. Whereas, the 
other two Syngnathidae presented diets with zooplanktonic 
affinity. Eel was distinguished by its benthophagous car¬ 
nivorous diet, while the European seabass presented pis¬ 
civorous affinity. The remaining species constituted a large 
multispecific cluster, which grouped the four Mugilidae spe¬ 
cies, as well as the barb, the big-scale sand smelt, the South 
European toothcarp and the sole. Mugilidae presented the 
highest inter-species similarity (Cbetween 0.8 and 0.98). 
The same observation was also noted for B. callensis and 
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Table II. - Abbreviations key of the several items (prey) identified in the stomach contents of the Ichkeul fish species. 


Prey groups 

Abbreviation 

Prey species 

Details 



Ab_alb 

Abra alba (probably Abra sp. and/or Scrobicularia plana) 

Bivalves 

Biv 

Cera_gl 

Cerastoderma glaucum 



Debrijbiv 

Bivalves debris ( Abra alba and/or Cerastoderma glaucum ) 



Hydrobia 

Hydrobia spp. ( H. ventrosa and/or H. acuta) 

Gastropods 

Gast 

Gast_debri 

Gastropods debris (probably Haminoea sp.) 



Moll ter 

Terrestrial molluscs (gastropods [snail]) 

Fishes 

Fish 

Fish 

Fish scraps (scales, fins, vertebrae, spines, etc.) 

Decapod 


Shirp 

Shrimp (probably Atyaephyra desmaresti incidentally juvenile of 
Parapenaeus longirostris and / or Penaeus sp.) 

crustaceans 


Crab 

Crabs [decapods] (debris, hooks, pincers) ( Carcinus aestuarii and / 
or Potamon algeriense) 



Pera_debri 

Peracarid crustaceans debris (amphipods and isopods) 



Ido_balth 

Idotea balthica 

Peracarid 


Sphaero 

Sphaeroma hookeri 

crustaceans 


Gamar 

Gammarus spp. (Gammarus aequicauda, Gammarus insensibilis 
and probably Gammarus locusta) 



Coroph 

Corophium volutator 

Worms 

Worm 

Annel 

Annelids ( Nereis diversicolor) 

Plath 

Plathelminths 

Insects 

Insect 

Insect 

Insects (insect scraps: hymenoptera wings [Formicidae], scraps, 
mandibles, coleopera hook [Blaps] or potentially Odonata) 

Mosqt 

Mosquito larvae [Culicinae] (larvae, nymphs and other forms) 



Chiron 

Chironoms (insects) 

Animal scraps 

Anm scrap 

Anim debri 

Various animal debris 



Cop 

Copepods (Oithonidae, Acartiidae, potentially Temoridae) 



Other_zooplk 

Other zooplankton (egg, larvae, veliger larvae, Artemia, nauplii etc.) 

Zooplankton 

Zooplk 

Mysid 

Mysid? (Unconfirmed, may be confused with a larval form 
euphausiids or other decapod) 



Clado 

Cladocerans (Zooplankton, Crustacea, Branchiopoda: several 
species) 

Zoobenthos 

Zoobent 

Ostrac 

Ostracods (Crustacea, zoobenthos) 

Foram 

Foraminifera (benthic protozoa) 



Rodo 

Rhodophyta ( Polysiphonia sp.) 



Rupp 

Ruppia maritima 

Macrophtyte 

Green algae 


Plant_seed 

Plant seeds (probably seeds of Potamogeton maritimus , incidentally 
seeds (akene) of Ruppia maritima or Scirpus maritimus) 

Mphyt/Alg 

Plant_debri 

Plant debris (Potamogeton maritimus and/or phragmites and 
potentially terrestrial plant debris) 



Alg_green 

Ulva spp. and Chaetomorpha spp., Cladophora sp.) 



Micro_benth 

([grazing]) microalgae and sediment superficial vegetal film 
(benthic grazing film ([grazing]) 

Other 

Unidentifed 

Other 

Other_[ Unident ] 

Unidentified; Organic aggregates. Mineral aggregats, calcareous 
annelid tubules annelids, others 

Sediment 

[Substratum] 

Sedim 

Sedim 

Sediment (mud / silt / clay / with sometimes a sand fraction in small 
proportions) 



nb FAS 

Number of full stomach analysed 



nb AS 

Number of analysed stomach 



nb ES 

Number of empty stomach 



Cv 

Vacuity coefficient 
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Figure 4. - Hierarchical Cluster Analy¬ 
sis (HCA) of fish species according to 
their diet. The threshold value of Horn- 
Morisita distance adopted to extract 
groups was 0.6. See figure 3 for species 
names. 



S. senegalensis (C/„„ = 0.91) as for A. boyeri and A.fasciatus 
(C hrn = 0.84) (Tab. I). 

Spatio-temporal diet variability 

The individuals caught in the West were distinguished 
by a less important vegetal fraction than those caught in the 
East of the lake. The eel food diet was marked by an inter- 
seasonal and an intersectoral variability (Fig. 5A). During 
the dry season, the HCA showed a clear inter-sector differ¬ 
ences; eel consumed more bivalves in the West than in the 
East part of the lake (Fig. 8). A.fasciatus showed a stable 
diet profile during the different seasons with equal contribu¬ 
tions of plant, zooplankton, zoobenthos, animal debris and 
amphipod/isopod crustaceans (Fig. 5B). 

The big-scale sand smelt A. boyeri exhibited low diet 
variability along the two seasons in the same lake sector. 
However, the contribution of micro-gastropods ( Hydrobia 
ventrosa ) was more important (> 20%) in the West, while in 
the East it is characterized by a more important amphipod/ 
isopod contribution ( > 50%) (Figs 5C, 8). For the barb, the 
wet season was marked by more important plant contribution 
than the dry season. The peracarid crustacean fraction varied 
between the West (16-18%) and the East (6-7%) (Fig. 5D). 

The bivalve contribution (12-14%) in the eastern sector 


(wet season) and the western sector (dry season), respectively, 
as well as that of insects (6-9%), in the eastern part (wet sea¬ 
son) and the western part (dry season), respectively, are to be 
underlined. In the eastern sector, B. callensis exhibited a food 
diet rich in macrophytes/algae, crustaceans and bivalves, while 
in the western sector the animal fraction was more important 
and more diverse with an affinity for gastropods. 

During the wet season, B. belone showed a diet essen¬ 
tially composed by insects and amphipods/isopods; the 
species was characterized by its affinity towards peracarid 
crustaceans and insects (Fig. 7). These two prey categories 
represented more than 80% of the species food diet. During 
the dry season, B. belone showed a diet with more impor¬ 
tant insect and fish contributions in the eastern than in the 
western part. D. labrax was characterized by an ichthyopha¬ 
gous diet. These prey proportions varied between 23% in the 
East sector-Dry season (ED) to 45% in the West sector-Wet 
season (WW). The animal fraction varied also according to 
season and sector sampled between 75% (ED) to 84% in the 
West sector-Dry season (WD) (Figs 5F, 7,8). 

A pronounced variability in the diet composition of 
S. senegalensis appeared; during the wet season, amphipods 
and isopods fraction was very important. During the dry sea¬ 
son, gastropods, plant and sediment fraction were dominant 
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Figure 5. - Spatiotemporal evolution of volumetric percentages (%V) of the different prey categories found in the stomach contents. 
A: Anguilla anguilla; B: Aphanius fasciatus; C: Atherina boyeri; D: Barbus callensis ; E: Belone belone ; F: Dicentrarchus labrax\ 
G: Engraulis encrasicolus', H: Solea senegalensis, in the Ichkeul Lake. ED: East sector-Dry season; EW: East sector-Wet season; WD: 
West sector-Dry season; WW: West sector-Wet season; n = N: number of non-empty analyzed stomachs for each species. 
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Figure 6. - Spatio-temporal evolution of volumetric percentages (%V) of the different prey categories found in the stomach contents. 
A: Liza aurata; B: Liza ramada; C: Liza saliens ; D: Mugil cephalus ; E: Pomatoschistus microps ; F: Syngnathus abaster; G: Syngnathus 
acus; H: Syngnathus typhle, in the Ichkeul Lake. ED: East sector-Dry season; EW: East sector-Wet season; WD: West sector-Dry season; 
WW: West sector-Wet season; n = N: number of non-empty analyzed stomachs for each species. 
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Figure 7. - Correspondence Analy¬ 
sis (CA) between fish species and food 
diets during the wet season. The addi¬ 
tive symbol attributed for each species 
referred to sectors; [E] for Eastern sec¬ 
tor and [W] for Western sector. 



Figure 8. - Correspondence 
Analysis (CA) between fish 
species and food diets dur¬ 
ing the dry season. The addi¬ 
tive symbol attributed for each 
species referred to sectors; [E] 
for Eastern sector and [W] for 
Western sector. 
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Figure 9. - Structural diagram of the Ichkeul Lake trophic web and trophic 
relationships (prey-predator) between ichthyofauna and their main prey 
(thick line = high predation; single line = ordinary to low predation; Dashed 
lines: larvae + juvenile of fishes (Predation fish to fish), thick dashed line = 
hight predation, thin dased line = ordinary to low predation). 


(Fig. 5H). The macrophytes and algae percentages 
varied between both lake sectors, but this variability 
was not the same for the four mugilid species. This 
prey variability was more important in the western 
part than the eastern one for L. aurata, when it was 
more important in the eastern sector than in the west¬ 
ern one for L. ramada and L. saliens (Fig. 6A-C). The 
diet composition was characterized by little changing 
percentages across sectors and seasons for M. cepha- 
lus (Fig. 6D). In the western sector, L. aurata seemed 
to ingest more detritus and sediment than in the east¬ 
ern one. L. saliens and M. cephalus were characterized 
by a diet with a low variability during the wet season 
for both sectors (Fig. 7). The diet of the Gobiidae P. 
microps was distinguished by the presence of six prey 
categories with same proportions for both seasons in 
the eastern part of the lake. 

The inter-seasonal variations of P. microps diet 
concerned the zoobenthos (between 18 and 38%), 
zooplankton (between 25 and 50%) and amphipod/ 
isopod crustaceans (between 16 and 22%) (Fig. 6E). 
The CA showed that diets of P. microps and .S', typhle 
were characterized by a zooplankton/zoobenthos 
affinity (Fig. 7). The two Syngnathidae S. abaster 
and S. acus presented a similar diet in both wet and 
dry seasons (Fig. 6F, G). Their diet was dominated by 
macrophytes/algae and zooplankton (Figs 7, 8). The 
spatio-temporal evolution of the different stomachal 
analysis indices ( nbAS , nbES, Cv ) of the different fish 
species is detailed in Appendix II. 

Figure 9 shows the main trophic relations linking 
the principal prey and fishes. It highlights the pre¬ 
dation relationships between fish species under the 
Ichkeul lacustrine ecosystem. The main prey were at a 
lower trophic level compared to fish. They were main¬ 
ly peracarid crutaceans, gastropods, bivalves, worms, 
zooplankton, zoobenthos, macrophytes and algae. 
Carnivorous species such as A. anguilla and D. labrax 
occupied particular positions and were characterized 
by the predation of both high (mainly fish) and low 
level preys, particularly small-sized prey and fry of 
other fish species with a lower trophic level (Appen¬ 
dix II). B. belone was characterized by a more specific 
and piscivorous diet (Appendix II). Syngnathidae, E. 
encrasicolus and P. microps occupied monospecific 
positions, while Mugilidae, S. senegalensis, B. callen- 
sis, A.fasciatus and A. boyeri represented multispe¬ 
cific groups in the fish food webs (Fig. 9). 
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DISCUSSION 

Environmental characteristics, sampling 
Environmental characteristics 

During this study, temperature and salinity values were 
typical to the studied environment and presented the same 
tendancy recorded in previous studies (ANPE, 2002-2007; 
Ben M’Barek and Slim-Shimi, 2002; Casagranda and Bou- 
douresque, 2010). In particular, salinity is characterized by 
important seasonal variations. Temperature tends to show 
spatio-seasonal variations and the eastern sector appears 
as the region presenting the lowest temperatures during the 
whole year, except during the dry months (July, August and 
September). This tend was already shown by Ben Rejeb- 
Jenhani (1989), Ben Garali et al. (2008) and Casagranda and 
Boudouresque (2010). This could be explained by the influ¬ 
ence of water currents introduced across the Tinja Oued and 
its role on marine-lagooner (Bizerta lagoon) water masses 
mixing in the eastern sector (Krai'em et al., 2003; Shili et 
al., 2007). The western sector and the central zone are less 
dynamic; they are shallower than the eastern area (BCEOM, 
1995; Ben M’Barek and Slim-Shimi, 2002; Ben Garali et 
al., 2008). The large macrophytal ( Potamogetonpectinatus, 
Ruppia cirrhosa) coverage contributes also to heat waters of 
this region, in particular during the dry session (but still less 
than eastern temperature during the dry period) when there 
is the decomposition of a huge quantity of these plants (Shili 
et al., 2007). 

The annual salinity oscillates between marine salinity 
during the hot season and freshwater salinity during the cold 
season (Ben Rejeb-Jenhani, 1989; Casagranda and Bou¬ 
douresque, 2010). The variability of the surface area and the 
water level are directly related to the water supplies routed 
to the lake and mainly to the region pluviometry (BCEOM, 
1995, Ben M’Barek, 2001; Ben M’Barek and Slim-Shimi, 
2002; ANPE, 2002-2007; Casagranda and Boudouresque, 
2010 ). 

Trophic diet 

The choice of volumetric analysis for point method did 
allow the analysis of fish diets when the microscopic organ¬ 
ism measurement (microalgae, zooplankton, etc.) is difficult 
and time-consuming (FAO, 1974; Hyslop, 1980). For many 
species, the low sample number (scarcity and difficulty to 
catch these species, particularly those of small size (Syn- 
gnathidae, P. microps, E. encrasicolus, A. boyeri ) allows 
only to get a superficial idea of their diets, and the results of 
their stomach content analysis should be taken with much 
caution. It is an important problem of this work. However, 
it is the first study which covers the whole trophic function¬ 
ing of the fish community in this lacustrine ecosystem. This 
study highlights the trophic relationships between prey and 
fish, and among fishes. This is an important step to charac¬ 


terize the food webs and understand the trophic functioning 
of this particular ecosystem. 

Prey like insects, micro-zoobenthos, percarids and benth¬ 
ic micro-gastropods are always present in the totality of the 
stomach contents analyzed, which characterizes the Ichkeul 
trophic network. In fact, the presence of these species is one 
of paralic ecosystem characteristics and were present in 
the same time in lacustrine, mudflats (swamp) and lagoon- 
marine habitats. 

From an abiotic point of view, these paralic environ¬ 
ments are characterized by shallowness, water stratification 
and high confinement that generate a high phytoplanktonic 
biomass (Frisoni et al., 1983; Guelorget et al., 1989) and 
high macrophytal biomass (Casagranda et al., 2005a; Casa¬ 
granda et al., 2006; Shili et al., 2007). The omnipresence of 
vegetal fraction within the stomach contents of numerous 
species attests the importance of the vegetal biomass in the 
lacustrine ecosystem functioning. 

This structural trophic network characteristics of Ichkeul 
Lake are similar to those observed in other similar environ¬ 
ments, like the neighbouring lakes and estuaries of North 
Africa (Guelorget et al., 1989; Bouchereau et al., 2000; 
Bazairi et al., 2005) or the lakes and transitional environ¬ 
ments of south Europe: Donana and Mar Menor (South of 
Spain) (Giordani et al., 2009), Ria Formosa (Gamito et al., 

2003) et Ria d’Aveiro (Garnerot et al., 2004) (South of Por¬ 
tugal), Fogliano Lake and Caprolace Lake (Tyrrhenian basin, 
central West, Italy) (Mariani 2001), Vaccares pond and Or 
pond (Camargue, Rhone Delta, South of France) (Mouil- 
lot et al., 2005 ; Vizzini et al., 2005 ; Dierking et al., 2012); 
and La Palme pond (Natural park of Narbonnaise, South of 
France ) (Carlier et al., 2008). 

The fact that a prey is found in the stomach contents of 
big-sized predator fish (A. anguilla, B. callensis , D. labrax ) 
and in those of small-sized ones ( B. belone, A. fasciatus, A. 
boyeri) suggests that food diet reflects more the availability 
and richness of prey in the environment than a direct influ¬ 
ence of fish size (morphological intrinsic parameter) or con¬ 
finement degree (physico-chemical abiotic parameter), as it 
was already observed in other Mediterranean lagoons such 
as Nador or Aveiro (Bouchereau et al., 2000; Garnerot et al., 

2004) . This availability was confirmed by the previous stud¬ 
ies in Ichkeul Lake, which pointed out the high biomass of 
hydrobies, percarid crustaceans and marcophytal biomass 
(Casagranda et al., 2006; Shili et al., 2007; Mancinelli, 
2012 ). 

According to Bazairi et al. (2005), the trophic organiza¬ 
tion of environments similar to Ichkeul is depending on a low 
number of species that dominate quantitively fauna commu¬ 
nities, particularly bivalves ( Cerastoderma spp., Abra spp., 
Scrobicularia plana) and amphipod/isopod crustaceans, 
when micro-grazers and herbivores are generally scarce, 
except some specific species such as hydrobies. The benth- 
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ic macrofauna is composed of strictly paralic species such 
as peracarid crustaceans (Corophium sp. Gammarus sp., 
Sphaeroma hookeri ), polychaetes {Nereis sp.), chironomid 
larvae and filter feeder pelecypods {Abra sp., Cerastoderma 
glaucum, Scrobiculariaplana), in addition to Hydrobia sp. 
This latter genus presents a quasi-monospecific population 
usually in the most confined zones (Guelorget et al., 1989; 
Lam Hoai and Amanieu, 1989; Garnerot et al., 2004). 

According to Frisoni et al. (1983) and to Guelorget et al. 
(1989), this benthic macrofauna is distributed all over those 
environments, depending on zones and confinement degrees. 
These are actually low diversified communities, with mixed 
strictly paralic species and some representatives of freshwa¬ 
ter fauna in the marginal remotes and outfalls of oueds fol¬ 
lowing the lake contour (Guelorget et al., 1989). Moreover, 
only a few crustaceans living in paralic environments con¬ 
stitute a principal prey for most of the Ichkeul fish species. 
Amphipods and isopods constitute an important biomass 
in this environment (Casagranda et al., 2006; Mancinelli, 
2008) and their high availability may explain their occur¬ 
rence in the food diet of fishes (particularly carnivores). 
Other groups can be also considered as main fish prey. Small 
meiofauna organisms (ostracods and foraminifera) (Guelor¬ 
get et al., 1999) are found regularly within digestive tracts of 
fish, from small sized strictly benthophagous species (Syn- 
gnathus spp., cyprinodon, etc.) to bigger piscivorous species 
(eel, seabass), as well as in detritivorous/omnivorous species 
(mullets and barb). These incidentally ingested preys are 
constantly present in almost all aquatic ecosystems (Zaibi et 
al., 2011). 

Spatio-seasonal variability 

The spatiotemporal variability of stomach contents was 
characterized for many fish species. This variability could 
be related to different factors: trophic preferences versus 
resource availability according to sector and season. The 
second factor is closely linked to the influence of salty and 
freshwater interface (marine water versus continental water) 
(Mariani, 2001). 

The structure of the lacustrine prey and fish community 
reflects the physical, geochemical and hydrodynamic al char¬ 
acteristics of the environment. The fish distribution is cor¬ 
related with marine influence degree on such lacustrine eco¬ 
systems as Ichkeul Lake (Guelorget and Perthuisot, 1983, 
1992; Mariani, 2001). Added to the constraints of these abi¬ 
otic parameters, one should also consider the opportunity for 
each fish species to find food adapted to its size and physi¬ 
ological needs. This relation between size, diet, prey and 
ecological conditions constitutes an important parameter to 
understand the space occupation by a species (Paugy and 
Leveque, 1999). 


Environment predator-prey relationship 

The important biomass of the prey categories (amphi¬ 
pods, isopods) and their predators (cyprinodon, sand smelt, 
anchovy, etc.) is due to the relative ‘stability’ of the varying 
ecological and trophic conditions, which are favourable to 
these species (water quality, physico-chemical parameters 
and nutrients availability). This can explain the low degree 
of spatiotemporal diet variability of predator fishes; the 
fact that Ichkeul Lake does not have a pronounced confine¬ 
ment gradient (ovoidal form, important freshwater supplies, 
strong marine influence and typical hydrodynamism) induc¬ 
es an important biomass of benthic macrofauna (Frisoni et 
al., 1983; Bouchereau et al., 2000; Garnerot et al., 2004). 

This hypothesis is supported by works of Mancinelli 
(2012) and Kuan-Yi et al. (2009) about the close relation¬ 
ship between macrophytal biomass and then - associated epi¬ 
phytes and epifauna, particularly the peracarid crustaceans, 
which were found amply into fish stomach contents. 

The bivalves were found particularly in stomach con¬ 
tents of eels during the dry season. These molluscs were 
found in important quantities in muddy or/and muddy-slimy 
lacustrine sediments (Gamito et al., 2003; Bouchereau et al., 
2009a, b). The benthophagous behaviour of eels was already 
demonstrated in Prevost pond, Ingril and Mauguio lagoons 
(Guelorget and Perthuisot, 1982; Bouchereau et al., 2009b). 

Small-sized fish species, such as Aphanius fasciatus 
and Atherina boyeri , represented an important part of the 
diet of piscivorous and carnivorous species ( D. labrax and 
A. anguilla). Figure 9 proposes a synthetic diagram of the 
trophic relationships evidenced in the present study that link 
fishes to their principal prey within the Ichkeul Lake ecosys¬ 
tem. 

Seasonality 

Seasonality of fish diet was not observed for all fish 
species sampled. A clear seasonality and a noticeable diet 
change, between dry and wet seasons, are particularly obvi¬ 
ous for A. anguilla, S. senegalensis and B. callensis. This last 
species is characterized by a high diet variability compared 
to the two other fish species, as it does not find the same prey 
in the freshwater environment of rivers, in which it ascends 
during the dry season, than those available during the wet 
season in the Ichkeul Lake. The availability of the lacustrine 
prey for the barb is especially concentrated in the western 
sector, which is more influenced by continental water sup¬ 
plies (Casagranda et al., 2005a, 2006). 

Thus, D. labrax, S. senegalensis, B. callensis and A. 
anguilla showed a spatio-temporal variability of their diet. 
The benthophagous diet of these species varies consider¬ 
ably between wet and dry seasons. In fact, the transition of 
a diet ( i.e zoobenthophagous to detritivorous) was already 
observed in shallow paralic environments, which are influ¬ 
enced by abiotic parameters, particularly salinity and tem- 
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perature (Paugy and Leveque, 1999; Casagranda et al., 
2005a, 2006; Casagranda and Boudouresque, 2010). 

Therefore, except these previously mentioned spe¬ 
cies, there is no trophic variability for the other fish species 
despite a high physico-chemical seasonality of the lacustrine 
ecosystem (salinity, temperature, water level, macrophyte 
coverage) (Casagranda et al., 2005a, 2006 Casagranda and 
Boudouresque, 2010). 

This can be explained by the abundance of resources 
during the whole year, which constitutes the main prey of 
most fish species. Some of those preys are constantly avail¬ 
able in the environment and their biomasses are very impor¬ 
tant. This obervation was confirmed by many authors like 
Guelorget et al. (1987, 1989) or Garnerot et al. (2004). 
Indeed, biomasses of zoobenthos prey species are high in 
paralic environments, while specific richness is low. These 
prey have a large tolerance and high ecological amplitude 
valence, which allow them to adapt to the physico-chemical 
variability of their environment (Mariani, 2001; Garnerot et 
al., 2004; Mancinelli, 2012). 

Species like mullets, syngnathidae, cyprinodons, that 
showed a low spatial and/or seasonal variability of their food 
diet, composed a distinct multispecific trophic group. These 
species presented more or less similar diet composition as 
they move easily from eastern to western lake sectors during 
wet and dry seasons. 

CONCLUSION 

Prey ingested by fish species in Ichkeul Lake belonged 
mainly to lacustrine benthic macrofauna, in particular gas¬ 
tropods ( Hydrobia spp.), bivalves (Abra alba, Cerastoderma 
glaucum) , and peracarid crustaceans (amphipods/isopods). 

The omnipresence of vegetal fraction into stomach con¬ 
tents of the majority of fish sampled demonstrates its trophic 
importance, as well as the importance of this biomass in the 
lacustrine ecosystem. The majority of fish trophic categories 
defined according to food diet are the detritivores and the 
benthophagous carnivores. 

The seasonal and spatial variability of food diet differed 
greatly from a species to another. Anguilla anguilla, Solea 
senegalensis and Barbus callensis presented the strongest 
variations. These variations were particularly related to the 
nature and the proportions of the ingested preys. 

The trophic clusters and the spatiotemporal food vari¬ 
ability of fish species reflected the lacustrine trophic network 
structure and its evolution. 

The hydrological changes, linked to the water resources 
rarefaction supplying Ichkeul Lake, leads us to question its 
repercussion on the possible structuring of the trophic net¬ 
works . 
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Appendix I. - Number of analyzed stomachs, number of empty stomachs, vacuity coefficient Cv (%) between sector and seasons for the 
several fish species. Cv: Vacuity coefficient; ED: East sector-Dry season; EW: East sector-Wet season; N: number of non-empty analyzed 
stomachs for each species; N t : number of analyzed stomachs for each specie; N e : number of empty analyzed stomachs; WD: West sector- 
Dry season; WW: West sector-Wet season. 
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25 

7 

28.0 

27 

13 

48.1 
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Appendix II. - Trophic level and classification by trophic category of the several species according to fishbase (www.fishbase.org) and 
Karachle and Stergiou, 2012. 


Species 

Trophic level * 
(www.fishbase .org) 

Trophic category** 
(Karachle & Stergiou, 2012) 

Trophic category (www.fishbase.org) 

Anguilla anguilla 

3.5 ±0.6 

OA/CD 

Carnivores mesophagous 

Aphanius fasciatus 

2.7 ±0.2 

OV 

Detritivores (Omnivores) 

Atherina boyeri 

2.3 ±0.3 

OV 

Carnivores microphagous 

Barbus callensis 

2.8 ±0.3 

OA 

Detritivores (Omnivores) 

Belone belone 

4.2 ±0.7 

CC 

Piscivores 

Dicentrarchus labrax 

3.8 ±0.6 

CC 

Carnivores mesophagous/macrophagous 

Engraulis encrasicolus 

3.1 ±0.4 

OA 

Carnivores mesophagous 

Liza aurata 

2.5 ±0.2 

OV 

Detritivores 

Liza ramada 

2.2 ±0.1 

OV 

Detritivores 

Liza saliens 

3.0 ±0.2 

OA 

Detritivores 

Mugil cephalus 

2.1 ±0.2 

OV 

Detritivores 

Solea senegalensis 

3.1 ±0.3 

OA 

Carnivores mesophagous 

Syngnathus abaster 

3.2 ±0.4 

OA 

Carnivores mesophagous 

Syngnathus acus 

3.4 ±0.6 

OA 

Carnivores mesophagous 

Syngnathus typhle 

4.3 ±0.8 

CC 

Piscivores 

Pomatoschistus microps 

3.3 ±0.4 

OA 

Carnivores mesophagous 


* Calculated from diet study species trophic category was identified using their trophic level t (www.fishbase.org) 
** Trophic category (Karachle & Stergiou, 2012): According to this classification; 


H = pure herbivores (2,0 <t <2.1); 

VO = omnivores with vegetal preference (2,1 <t <2,9); 

OA = omnivores with animal preference (2,9 <t <3,7); 

CD = carnivores with preference for decapods/fish (3,7 <t <4,0); 
CC = carnivores with preference for fish/cephalopods (4.0 <t <4,5). 
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